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Abstract 

We explore the feasibility of using high conversion-efficiency periodically-poled crystals to pro- 
duce photon pairs for photon-counting detector calibrations at 1550 nm. The goal is the develop- 
ment of an appropriate parametric down-conversion (PDC) source at telecom wavelengths meeting 
the requirements of high-efficiency pair production and collection in single spectral and spatial 
modes (single- mode fibers). We propose a protocol to optimize the photon collection, noise lev- 
els and the uncertainty evaluation. This study ties together the results of our efforts to model 
the single-mode heralding efficiency of a two-photon PDC source and to estimate the heralding 
uncertainty of such a source. 

PACS numbers: 
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I. INTRODUCTION 



Parametric down-conversion (PDC) consumes pump photons and produces light with a 
two-photon field description [1]. This two-photon light, which allows one photon to indicate 
or herald the existence of its twin, has a key role in applications such as quantum radiometry 
Q], where it makes possible an independent primary standard method for direct photon- 
counting detector calibrations. In addition, this two-photon light is key to improve single- 
photon source schemes that have significant advantages over the attenuated laser sources 
often employed in quantum information applications. To realize a single-photon source 

appropriate for these applications, we consider two key issues- a) efficient photon pair 
production and b) collection of the heralded photon in a single- mode fiber, which is of interest 
for telecom applications. Periodically poled crystals (PPLs) have been used for efficient pair 
production p, but have not yet been used for metrology. While metrology applications 
would certainly benefit from the high pump-conversion efficiency of PPLs, characteristics 
such as the level of background fluorescence may limit the uncertainty that can be achieved. 

As for the collection efficiency of PDC light into a single-mode fiber, the goal has been 
to understand how the detection of one photon of a PDC pair in a well defined spatial 
and spectral mode, defines its partner's spatial and spectral mode that can then be effi- 
ciently collected. Recent efforts have included experimental and theoretical work in both 
continuous wave (CW) 0,0,0] and pulse-pumped 1(3, 11, Q| bulk crystal PDC, and have 
highlighted the critical problem of coupling the PDC photons into single-mode fibers (SMF) 
due to the difficulty of spatial and spectral matching of the photon pairs. The results show 
that optimization of the heralding efficiency is very sensitive to spectral and spatial mode 
selection. 

Here we describe a PDC source and a measurement technique and compare it to other 
configurations. We present theoretical and experimental characterizations of the photon 
collection problem and suggest possible solutions. 

II. EXPERIMENTAL SETUP 

We pumped (Fig. |TJ) a 5 mm long periodically poled MgO-doped lithium niobate (PPLN) 
crystal with a CW laser at 532 nm. We used noncritical phase- matching (90° phase-matching 
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angle) with a 7.36 /iin poling period to produce 810 nm and 1550 nm photon pairs at external 
angles of 1° and 2°. We achieved fine tuning by adjusting the crystal temperature near 131° 
C. Lens L p in the pump beam produced a gaussian beam waist of w p ~144yum at the crystal; 
cutoff filter Fc blocked the pump laser; and dichroic mirror DM separated the 810 nm (beam 
1) and 1550 nm (beam 2) photons. An extra interference filter (F x ) at 810 nm with a full 
width half-maximum (FWHM) of 10 nm suppressed extra fluorescence from the PPLN to 
reduce background heralding counts, while the SMF geometry imposed spectral selection 
with a heralding bandwidth, Ai of ~2 nm FWHM (see Table I). L 12 were aspheric coupling 
lenses of focal length 8 mm, AR-coated for 810 nm and 1550 nm. A monochromator in the 
heralding arm measured its spectral width. The heralding arm was routed to a SMF and then 
to Di, a Si Avalanche Photodiode (APD), while the heralded arm also coupled to a SMF, 
was sent to D2, an InGaAs APD operating in gated mode. We gated the InGaAs detector 
with the detection of an 810 nm photon in the heralding arm and adjusted the delay of the 
heralded arm with appropriate length SMF and electronic cables. Calibrating the InGaAs 
detector, including the SMF and coupling lens optical losses against a conventional detector 
standard using an attenuated laser source, yielded a detection efficiency of rj^et = 9.8%±0.5% 
(k=l, absolute uncertainty). 



III. MEASUREMENT TECHNIQUE 

The heralded D 2 detector efficiency as measured by the PDC source is given by 

Vdet = W (XP • r opt • ^SMF-lcns) (1) 

where xp is the heralding efficiency (single-mode preparation efficiency, discussed in the next 
section) j^Q], r op t is the overall heralded arm optical transmittance (including PPLN, Fc, 
and DM), and rsMF-iens is the optical transmittance of L 2 and the SMF on the heralded arm. 
Xd is the raw D 2 channel detection efficiency, directly measured according to 

P — P 
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where P CO inc is the probability of coincidence counts per gate, P unC orr is the probability of 
uncorrelated or accidental coincidence counts per gate (determined by changing the heralding 
delay so the detection gate misses heralded photons), and ^backgnd g i s the probability of 
gating counts produced by uncorrelated photons and dark counts in the heralding arm. 



A. Raw detection efficiency 



To understand Eq.(J2J) we note that when Di fires, D2 is gated for a duration T. We define 
the probability of D 2 not firing during T as pj.. We assume that the "events" that make D 2 
fire have a poissonian distribution (ie., uniform in time). D± can fire for either background 
or PDC photons, and the associated probabilities are P^^^ g and -P P dc ms = 1 — P D ackgnd S - 
The probability that D 2 fires for true correlated photons is -P C o! n c ldcd = PpDc dmg PT/2XD, where 
Pt/2 i s ^he probability for the heralded detector to fire for the correlated photon arriving 
at T/2, as opposed to firing during the first half of the gating time. (We assume that 
correlated photons arrive exactly at T/2.) The probability that the heralded detector fires 
for accidental events is 

Puncorr = (1 - + C 1 " ^^PDC^ + Pt/ 2 (1 ~ Xd)(1 - Pt/^PDC^ (3) 

The first term is the probability of an accidental coincidence due to Dx firing for a background 
event and D 2 firing for any event. The second and third terms are for Di firing for a PDC 
event and D 2 firing due to a background event either before or after the the arrival of the 
PDC photon at T/2. For D 2 to fire after T/2 it cannot have fired due to a background 
photon before T/2 or a PDC photon at T/2. Those factors make up the third term. The 
total coincidence probability is P CO inc = P C oinc lded + -Pmcon- Considering that for poisson 
statistics pj, = p\-i 2 x Pt/2' anc ^ inverting the above equations, we obtain Eq.(J2J), where 

-Puncorr 1 Pt~ 



B. Estimate and uncertainty 



We develop an estimate of the raw detection efficiency and the associated uncertainty 
following Ref . . The terms P C oinc, -Puncorr, and P^l^^ g in Eq.(2) are all independent 
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statistical variables, therefore the estimate of the raw detection efficiency is 
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Applying the maximum likelihood model estimator to the random variables P C oinc and -P U ncom 
the probability of M coinc coincidence counts given Mh er aiding heralding counts is 

P(M coinc |M heraldmg ,p) = - }U f he ™ ldm i M r. x P M — (1 - p)^, s ~M coin% (5) 

-'"coinc • ^"heralding -'"coinc J- 



where p is the parameter to estimate. The maximum likelihood function is then 

L{M co i nc , -^heralding \p) = -P(^coinc | -^heralding; p) x P ( ^heralding ) j 

where -P( Mheraiding) is the distribution probability of heralding counts, which is not known. 
The minimum of the function £(M coinc |M heralding ,p) gives the estimate of P C oinc, without the 
need of -P(M ne raiding)- It can be shown that (P C oinc) = (M co i nc /Mh era i drng ). Less trivial is 
the estimate of (P U ncorr)- We must estimate 1/(1 — P U ncorr), as is clear from Eq.(jU), so we 
introduce the parameter p' as p = (p' — l)/p'- The maximum likelihood function is in this 
case 

L ( M uncorr , M hera f ding \p ) = P ( M uncorr | M hera f ding ,p)P( ^herafding ) ) ( 6 ) 

which we minimize to obtain ((1/(1 - P U ncorr)) = (A^SdJg/ '{ M heSn g ~ Mmcorr))- We 
distinguished M heralding from M her *^ g because we measure heralding counts in two configu- 
rations, one with the presence of true coincidence photons and one in absence of coincidence 
when an electronic delay line is applied to the heralding arm. M uncorr are the coincidences 
corresponding to uncorrelated events. For ^ackgnd 8 ^ ne maximum likelihood model is not 
applicable because we cannot measure it at the same time as the background counts M" backgnd 
and M hera i ding , therefore we apply (P^l^ g ) ~ (M backgnd )/(M hcra i ding ). This implies that 
(^backgnd^) mus t be close to 0, to have no impact on the estimate of the raw detection 
efficiency We finally estimate 
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heralding 

with the uncertainty given by gaussian uncertainty propagation for each random variable 
considered independently. 



C. Heralding efficiency 



The heralding efficiency xp is the efficiency of preparing a photon in the heralded channel 
in a definite spectral and spatial mode, by specific mode selection of the heralding arm. It 
quantifies how well the collection system geometrically catches photons correlated to those 
seen by Di. To calibrate a SMF-coupled detector, the heralding efficiency must be optimized 
and estimated. It has been estimated in ref. ^ for a bulk crystal. For PPLN k| it is 
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where io ,i,2 are the mode waists of the fibers at the crystal, and Ai^ are the FWHM of 
the spectral distribution selected geometrically by the SMFs. / is a correction function 
accounting for crystal length (L) and PDC crystal frequency dispersion given by [3| 



dx e - px2+3 ^ (Erf[^l - Erff^^l) 
fip, q) = ~ Lj ^ L A 0) 

The parameters for this case are: 

L 2 (w 2 a 2 s + wl 2 tan6> 2 + w^(a s + tanft) 2 ) 
C\ = C2 H f — ^ — '-7. 9 — 7 — 5- — 9 — r; ; (10) 

[<2™p + <l(™p + < 2 )] 
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Here, a = 2a/1u(2) converts between the FWHM and the gaussian profile 1/e 2 radius. The 
other terms are D\ s = .Di(cos#i — sin 6\ tan 9{) — D s (cos6 s + sin 6 S tan 6 S ), D P i = —Di(cos6i + 
sin 6i tan ft) + A>, with A, s = ^ffi^ ln,,., ^ P = d ^;^/ c [ np . a s = -cos# s tanft + 
shift,. ft )S are the associated idler and signal emission angles around <p = it/ 2 (the phase- 
matching angle in a non-critical phase- matching configuration), and Ti- hSjP (ui )S)P , O ) are the 
indices of refraction (all are e-rays) at the three frequencies. fii, s , p are the central angular 
frequencies. 

For our experiment the fiber mode field diameters are MFD!=3.9 /im for the heralding 
arm and MFD 2 =5.6 /xm for the heralded arm, giving mode waists u? ,i,2 = ^1,2 MFD 12 at the 
crystal according to the magnification (Mi^) used. Because of this, we identify a SMF as a 

spectral filter with a gaussian spectral distribution given by J S) i(o; Si i— Sj i) oc e Al 2 
The component of spectral width due to geometric selection, A 12 is given by the FWHM 
angular collection Aft j2 = a ^ 1 '" - and the spectral/angular spread of the PDC, ft iS (A Sj i) 
around the central wavelength A; iS . Here, because the non-degenerate PPLN configuration 
and the bandwidth estimation are critical, we measured it with a monochromator in the 
heralding arm. We spectrally scanned the heralding single counts at three different values of 

^ (TaMe °' We al8 ° theoretica " y ffi ™rf 86 bandwldths by evaluatmg the FWHM 

of the phase-matching function according to 15], extended to the case of PPLN. Here we 
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consider the PDC phase-matching function including the pump transverse distribution 

*K ft, ft) = exp p - y - x — — — . (11) 



V Ak z L 

A&; X) y jZ (u; s , ft) ft) are the mismatch of the k-wavevectors A/c x>yiZ = (k p — k s — ki) x y z , in terms 
of the pump, signal, and idler k-vectors 

n(u; p Vp n(o; B )a; s ra(a;i)a>i 27r 

Afc z = cos ft cos ft — 

c c c A 

AA; xy = ^^smft + ^^sinft. 
c c 

A is the poling period. To estimate the bandwidth selected by the heralding arm fiber, 
we first find Of**** that maximizes Eq. (|ll|) . and then we consider for each wavelength the 
function f{9 s ,u s ) = 1 for $(^ s , ft max , ft) > 0.5 and /(0 B ,w B ) = for $(w B , ft max , ft) < 0.5. 
Finally we invert the functions /(ft + A9x/2,u s ) and /(0 S — A6i/2,u s ) to determine Ai. 

We estimated the heralded arm bandwidth in the same way. This estimation, which 
was found to be in agreement with the experimental values, allowed us to extrapolate the 
bandwidths of the heralded arm for other configurations. The discrepancy between the 



estimated and measured values is likely due to using Sellmeier's equations for undoped 
PPLN. Moreover, we do not have precise knowledge of the crystal poling region length (it 
maybe slightly shorter than full 5 mm crystal length), a crucial engineering parameter. In 
Fig. |21 we report the estimated Xp versus w 0) i,2 f° r our experiment's fixed pump waist, 
and crystal length, in two spectral-selection configurations on the heralding arm, (a) the 
SMF acts as a spectral filter with an average FWHM of about 2 nm and (b) a FWHM 
0.1 nm spectral filter (a monochromator was used for this narrow bandwidth). The two 
experimental conditions for the measurements in Table II are roughly indicated on Fig. Efa) 
by the squares, while the conditions for the narrow band measurement is indicated on Fig. 
I2jb) . While the heralded arm is spectrally selected only by the fiber, the large improvement 
achieved by narrowing the heralding bandwidth indicates that spectral mode matching is 
more critical than spatial mode matching (the only spatial requirement is that w o 2 = w p , 
w ,i > w p ). Using this we should be able to improve Xp by further narrowing the heralding 
arm bandwidth. 
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IV. EXPERIMENTAL RESULTS 



Following the approach in section III.B, we measured the raw detection efficiency in our 
setup at the best signal-to-noise ratio of the heralding arm. Table II shows xt> measured with 
u> Q) i =82 /jm fixed (g?i=270 mm) for 2 different positions of the heralded arm, corresponding 
to w D) 2 =158 (197) (g?2=300 (380) mm). The estimates and uncertainties calculated 
with this model and the measured intrinsic statistical fluctuations show good agreement. 
For each setup we tested the repeatability by applying our alignment procedure Q over 
several days. Fig. 3 shows the measurement repeatability in agreement with the estimated 
uncertainty at the 2 o"ml level. However, we point out that the repeatability test was 
limited to collection lens positions close to the source. We obtained theoretical values of 
X P = 20%, and 13%, respectively for the waists in Table II using only the SMF spectral 
selectivity. The decrease of Xp for bigger u> 0i2 is in agreement with the measured decrease 
of Xd- Our theoretical model for Xp is very sensitive to the measurements of the beam 
waists and the bandwidth estimates. Estimating the optical transmittances at T opt = 65% 
and TsMF-icns=83% by an independent calibration and using the measured xd in Table II, 
we were able to consistently obtain Xp — 48%, 37%. On occasion we measured Xp much 
closer to one by inserting a monochromator in the heralding arm, however the stability of 
our setup did not allow repeatable results at this level. 

At this stage, we note that the model can only be considered qualitative due to the spectral 
approximation adopted. In our approach we limited the calculation of the mismatch term 
to first order in the transverse k-vector and the frequency, losing some of the correlation 
between the frequency and the spatial variables in the biphoton field Q|. To compensate, 
we introduced a selection term due to the SMF spatial selection, but that works best for 
the degenerate case where the PDC emission is broader and this correction has less of an 
effect on the final mode matching than in the nondegenerate case. That case requires a full 
numerical solution. 

V. CONCLUSIONS 

We have shown the feasibility of using PPLN for single-photon detector calibration at 
1550 nm, as well as for a single-photon single-mode source. We highlighted the need to 
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maintain low dark counts in the heralding arm to reduce statistical fluctuations. We have 
shown that proper spatial mode matching and spectral mode selection in the heralding arm 
is of paramount importance in achieving the highest heralding efficiency. 
This work was supported in part by ARDA, ARO, and DARPA/QUIST. 
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VI. FIGURES AND TABLES 



TABLE I: Measured and estimated Ai 
di w Q> i Ai Ai 

measured estimated 



(mm) (/mi) 


(nm) 


(nm) 


270 82 


2.05 


2.6 


400 105 


1.87 


2.3 


520 132 


1.64 


2.1 
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TABLE II: Measured xd with the estimate, maximum likelihood (ML) uncertainty and statistical 
fluctuations 



^0,2 =158 /im 


w ,2 =197 /im 


XD CML Cstat 


XD 0"ML Cstat 


(%) (%) (%) 


(%) (%) (%) 


2.508 0.044 0.040 


2.015 0.030 0.028 


2.555 0.025 0.019 


2.056 0.028 0.035 


2.584 0.036 0.031 


2.107 0.034 0.028 


2.489 0.039 0.037 


1.996 0.032 0.026 
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FIG. 1: Setup to herald single-photons from CW PDC from a PPLN crystal. 
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FIG. 2: Calculated \p versus heralding and heralded waists (a) without an interference filter 
and (b) with spectral filtering on the heralding arm of 0.1 nm FWHM. Estimated values in the 
experimental conditions are given by the black squares. w p = 144/zm for both figures. 
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FIG. 3: Repeated measurements of xd for 2 setups (k=2 ML uncertainty shown). 
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